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Abstract 

The clusters of up spins of a two-dimensional Ising ferromagnet undergo a second order 
percolative transition at temperatures above the Curie point. We show that in the scaling 
limit the percolation threshold is described by an integrable field theory and identify the non- 
perturbative mechanism which allows the percolative transition in absence of thermodynamic 
singularities. The analysis is extended to the Kertesz line along which the Coniglio-Klcin 
droplets percolate in a positive magnetic field. 



1 Introduction 



The Ising model provides the fundamental example of critical behavior produced by the interac- 
tion among an infinite number of degrees of freedom, represented by spins on a lattice. Another 
type of criticality, which requires no interaction and, for this reason, is often called geometri- 
cal, finds its simplest illustration in random percolation [1]. Here the sites of the lattice are 
randomly occupied with probability p, and criticality (the appearance of an infinite cluster of 
nearest-neighbor occupied sites) arises when p reaches a lattice-dependent value p c . 

On the other hand, a ferromagnetic system like the Ising model encodes itself a percolation 
problem: nearest-neighbor sites with spin 'up' form clusters which can percolate for specific 
values of the parameters (temperature T and magnetic field H) which determine the interaction 
among the spins. This kind of percolation problem, which is completely determined by the 
magnetic properties and does not affect them, is called correlated percolation. The comparison 
of the magnetic and percolative phase diagrams on the T-H plane is particularly interesting. It 
turns out that in two dimensions (not in three) the temperature T p at which clusters percolate 
in zero field coincides with the Curie temperature T c Even in this case, however, a second 
order percolation line goes from (T, H) = (T c , 0) to a finite value of the magnetic field at infinite 
temperature, a circumstance which can appear paradoxical if one considers that there are no 
thermodynamic singularities above T c and that spin-spin correlations decay exponentially. 

Essential insight into this problem comes from the Kasteleyn-Fortuin representation [3] which 
associates the percolative properties to an auxiliary site variable Sj which takes q values. As a 
result, the magnetic and percolative transitions, being related to different site variables, are in 
a sense disentangled. The price to pay is that Sj, being auxiliary, cannot remain in the game to 
the end: everything needs to be evaluated in the limit q — > 1 in which the Ising spins are the 
only real degrees of freedom. The solution of the above paradox then requires to deal with the 
subtleties of this limit. 

Field theory, as the natural framework for dealing with second order critical points and the 
scaling region around them, is in principle the right place where to address this problem and 
isolate its universal features. The difficulty, however, is in the fact that the answer can hardly 
come from a perturbative approach. In this paper we show that the percolation line above T c 
in two dimensions actually corresponds, in the scaling limit, to an integrable field theory. On 
its exact solution, we can perform the limit q — ► 1 analytically, unveiling the presence of q — 1 
massless particles together with a massive unstable particle with lifetime inversely proportional 
to q — 1. Percolative properties are determined at first order in q — 1, where the theory is 
massless and describes the crossover from the correlated to the random percolation fixed point. 
The magnetic properties are determined instead at q = 1, where there are no massless excitations 
left and the massive particle has become stable, providing the required finite correlation length. 

Most of this discussion for the clusters can be repeated for the Coniglio-Klein 'droplets' [lj 
which, for H = (also in three dimensions), satisfy the requirement of Fisher's droplet model 
[5]: T p = T c and percolative exponents equal to the magnetic ones. Also the Coniglio-Klein 
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droplets exhibit for H > a second order percolation line, known as the Kertesz line [6], going 
this time from the Curie point to a finite value of the temperature at infinite magnetic field. 
In the scaling limit this line corresponds again to a renormalization group trajectory within 
the Ising field theory, which, however, in this case is not integrable, so that its origin through 
the resonance mechanism observed for the clusters, although very likely, cannot be followed 
analytically. 

The paper is organized as follows. In the next section we review the main results about the 
phase diagram of Ising clusters before recalling in section 3 the Kasteleyn-Fortuin formulation 
of the problem and the renormalization group analysis which leads to the identification of two 
different fixed points for clusters and droplets. The field theoretical discussion of clusters and 
droplets is then presented in sections 4 and 5, respectively, while section 6 contains a summary 
of the main conclusions. 

2 Percolation of Ising clusters 

Consider the ferromagnetic Ising model defined by the reduced Hamiltonian 



where a% is a spin variable located at the i-th site of an infinite regular lattice, T > and H 
are couplings that we call temperature and magnetic field, respectively, and the first sum is 
restricted to nearest-neighbor spins. For H = the model is well known to exhibit a non-zero 
magnetization per site M = (<7j) at temperatures below a critical value T c (the magnetization 
in zero field is called spontaneous). The magnetization has a discontinuity at H = along a 
path taken at fixed T < T c on the T-ii-plane (first order transition), and vanishes when T c is 
approached from below at H = (second order transition). No other magnetic transition (i.e. 
discontinuity in M or its derivatives) takes place away from T < T c , H = 0. 

If, given a spin configuration, we draw a link between nearest-neighbor 'up' spins (i.e. spins 
taking the value +1), we obtain connected sets of up spins that we call clusters. Of course, 
clusters of 'down' (—1) spins are defined analogously; in the following, talking of clusters without 
further specification we will refer to clusters of up spins, being understood that similar statements 
hold for the clusters of down spins under the substitution H — > —H. 

For H = +oo all the spins are forced to be up, so that the whole lattice is occupied by a 
unique infinite cluster. The fraction P of the lattice occupied by this infinite cluster decreases 
when the magnetic field decreases at fixed T, and is certainly zero at H = — oo. We denote by 
Hq{T) the value of the magnetic field below which the infinite cluster is absent. Since at T = 
all the spins are up for H = + , we have Hq(0) = + . When T = oo, on the other hand, the 
spins are uncorrelated and take the value +1 with probability e^/2 cosh i7 ; hence, denoting by 
p® the percolation threshold for the random site percolation problem, we have 
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Like T c , p® is non-universal, i.e. depends on the structure of the lattice. Notice for example 
that, if p® < 1/2, H (oo) is negative, so that at sufficiently high temperature an infinite cluster 
of up spins exists also in a negative field and coexists with an infinite cluster of down spins in 
an interval of H around zero. 

For finite temperatures the Ising model encodes a generalized percolation problem in which 
the sites are not independent but interact through the Hamiltonian ([1]). The probability that a 
site has spin up is (M + l)/2 in terms of the average magnetization per site. Some early studies 
on Ising clusters can be found in [101 [TTl [21 [12] . Here we summarize the main evidences 
relying on the following three statements: 

i) Hq(T) is a monotonic function; 

ii) the existence of a spontaneous magnetization implies the presence of an infinite cluster; 

iii) let p be the probability that a site has spin up at H = 0, and p c the critical value above 
which an infinite cluster appears. Then p c does not exceed the value p® of random percolation. 

Statement i) can be justified observing that, as a consequence of the ferromagnetic inter- 
action, the fraction of up (down) spins increases as we decrease the temperature for a fixed 
positive (negative) value of the magnetic field, and eventually becomes 1 at T = 0. Hence, it is 
reasonable to expect that, if Hq(oo) is positive (negative), the strenght of the positive (negative) 
field needed to produce (destroy) the infinite cluster of up spins decreases with the temperature, 
until it vanishes (it is zero at T = 0). 

Statement ii) is a rigorous result of [2]. Notice that the opposite is not true: for Hq(oo) < 
there is an infinite cluster at large enough temperatures in zero field, but there is no spontaneous 
magnetization above T c . Also, magnetization at H ^ does not imply an infinite cluster: if 
< H < Hq(oo), there is no infinite cluster at infinite temperature in spite of the positive 
magnetization. 

Statement iii) is an observation of [10]. A way of understanding it is to think of the ferro- 
magnetic interaction as inducing an attraction among the clusters, which favors (with respect to 
the random case) the formation of larger clusters, and eventually of the infinite cluster. Indeed, 
if Zj s i n g is the partition function, the probability of a configuration is e^^^Q /Zj s i ng . While 
configurations with a fixed number of up spins are all equally probable in the random case, in 
presence of the interaction cluster formation lowers the energy and increases the probability. 

We can now distinguish three cases according to the value of the random percolation thresh- 
old p° c : 

a) p® > 1/2, i.e. Hq(oo) > 0. In this case i) prevents the coexistence of infinite clusters of up 
spins with infinite clusters of down spins. Hq(T) takes the value + at T = and is forced by 
ii) to stay constant up to T c ; above T c it is free to increase and to reach its positive asymptotic 
value (Fig. 1). So T c is also the percolation point in zero field; since the magnetization vanishes 
at this point, one has p c = 1/2, which is consistent with iii). 

b) p® < 1/2, i.e. Hq(oo) < 0. In this case Hq(T) is zero up to some value T p of the 
temperature, above which it becomes negative, so that an infinite cluster of up spins coexists 
with an infinite cluster of down spins in the region \H\ < —Hq(T). At H = 0~ , the fraction 
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Figure 1: Qualitative phase diagram expected when the critical probability p c of random site 
percolation is larger than 1/2. The thick line is the curve Hq{T) above which there is an infinite 
cluster of up spins. Its dashed portion indicates that the percolative transition is first order 
below T c . The Kertesz line (dotted) will be discussed in section 5. 

p of sites with spin up, which is zero at T = 0, increases with the temperature and reaches 
the critical value p c for T = T p . Since p = 1/2 for T > T c , iii) requires T p < T c . Hence, for 
T p < T < T c there are infinite clusters, with different origin and different density, on both sides 
of H = 0. 

c) p° = 1/2, i.e. Hq(oo) = 0. In this case Hq(T) = at all temperatures. Continuity with 
the previous two cases requires that T p — > T~ as Hq{oo) — > 0~, so that, like in case a), the 
percolative transition is first order below T c and, presumably, second order above. 

This classification according to p® leads to a distinction between the two- and three-dimensional 
cases [10] . Indeed, the two-dimensional lattices hava^ pj? > 1/2 |13|. 114]. and then fall into the 
cases a), c), which share the same critical pattern with both the magnetic and the percolative 
transition taking place at T c in zero field. For the square lattice Ising model this was confirmed 
by series expansions in [12] and proved rigorously in [2|. For three-dimensional lattices the ev- 
idence is that p® < 1/2 [13 [16], so that they fall into the case b), for which the magnetic and 
percolative transitions in zero field are not simultaneous. For the simple cubic lattice this was 
first seen in [8] using Monte Carlo simulations. 

3 Kasteleyn-Fortuin representation 
3.1 Dilute Potts model 

The correlated percolation problem introduced in the previous section admits a formulation in 
terms of auxiliary Potts variables which generalizes that originally given for random percolation 
by Kasteleyn and Fortuin [3]. 

Let us first of all rewrite, up to an inessential additive constant, the Hamiltonian ([1]) in the 

'in a given dimension d, the critical probability p° of random site percolation decreases as the coordination 
number C of the lattice increases. In d = 2, the triangular lattice (C = 6) has p° c = 1/2; in d = 3, the diamond 
lattice (C = 4) has p° ~ 0.43 (see pQ ). 



4 



o o 



• — • o o • 
o o • — • o 
•— — • — • o • 
• — 1 o • — • 

Figure 2: Site-bond configuration on the square lattice. The full (empty) circles correspond 
to occupied (empty) sites of the lattice gas. Bonds (thick segments) between nearest-neighbor 
occupied sites are present with probability ps and define the KF clusters. For pb = 1 all 
segments are thick and the KF clusters coincide wiht the Ising clusters. 



lattice gas language 

- Hising = 7^ E lit i + A E ^ ' U = (<ri + l)/2 = 0,1, (3) 

(ij) i 

where A = 2H — a/T, with a a lattice-dependent constant. Then consider the dilute g-state 
Potts model with Hamiltonian 

-H g = -Hising + J ^ tjtj (S SiiS . - 1) , Si = 1,2, . . . ,q. (4) 

(ij) 

The partition function allows the following representation in which the sum over the Potts 
variables Sj is replaced by one over bond variables [17] 

= EE^ 

{U} {Si} 

= Y. e ~ Hl8in9< l Nv Y.P h B^-PB) l q Nc , PB^l-e~ J , (5) 

{ti} G 

where N v is the number of empty (ti = 0) sites, the last sum is performed over the graphs 
G obtained putting bonds in all possible ways between nearest-neighbor sites belonging to the 
restricted lattice formed by the occupied (ti = 1) sites only, b is the number of bonds in the 
graph G, b is the number of absent bonds on the restricted lattice, and N c is the number of 
connected components in G. Such connected components are called Kasteleyn-Fortuin (KF) 
clusters (Fig. 2). 

For q = 1 there are no Potts degrees of freedom (Z\ = Zj s i ng ) and the second sum in © 
gives 1, showing that pb is simply the probability that a bond is present on the restricted lattice. 
By construction, the KF clusters live on the clusters of up spins of the Ising model, and become 
the Ising clusters for pb = 1. 

The number q of Potts colors appears in ([5]) as a parameter (which can be taken continuous) 
through power-like terms accounting for the fact that both the empty sites and the KF clusters 
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can take q colors. If X is an observable, the average associated to the partition function ([5]) 
with q = 1 is 

W = Z7L 9 E e ~ Hlama E X Pb( 1 - PB? ■ (6) 

{U} G 

When H = +oo there are no vacancies and the average reduces to the weighted sum over bond 
configurations on the whole lattice, i.e. the usual random bond percolation problem. For generic 
values of T and H, instead, © corresponds to a generalized percolation problem in precence 
of vacancies whose distribution is weighted by the Ising lattice gas Hamiltonian. For example, 
for X = N c , © gives the mean cluster number in this generalized percolation problem, as a 
function of pb- Evaluated at ps = 1, this quantity is the mean number of spin up clusters in 
the Ising model at the given values of T and H. 

Notice that, while the correlations among Ising spins affect the associated bond percolation 
problem, the opposite is not true. At q = 1 the Hamiltonian (J3J reduces to TCising and the 
parameter J (or ps) plays no dynamical role: the conjugated bond variables over which the 
second sum in © is performed only serve enumeration purposes. In particular, an observable 
X which does not depend on the bond variables goes out of the second sum, which then gives 
1, leaving us with the usual Ising thermodynamic average. 

The role of the parameter q is further clarified if we consider the Hamiltonian TC q obtained 
adding to TC q the term H X^*i(^st,i ~~ 1)> where the magnetic field H conjugated to the Potts 
variable is usually called 'ghost' field. The new partition function Z q is then given [T7] by Z q 
with q Nc replaced by Y[ r (<7 — l)e HSr + 1 , where r labels the KF clusters in G and S r is the 
number of sites in the r-th cluster. When we expand the free energy per site around q = 1, 

f q = -l\nZ q = Using ~ (q - l)F + 0((q - l) 2 ) , (7) 

the function dF/dH depends on H through terms containing the factor \/N ^ r S r e HSr = 
n sSe HS , where ns is the number of clusters of size S per site, so that 

i.e. F is the generating function for the moments of the cluster size distribution. The probability 
that a set of sites belongs to the same cluster can be otained along the same lines introducing a 
site-dependent ghost field [17]. The important message of (J7J) and (JSj) is that the dilute g-state 
Potts model coincides with the Ising model at q = 1 and determines the properties of KF clusters 
within the lattice gas at first order in q— 1; it describes Ising clusters at first order in q — 1 when 
Pb = 1- 

3.2 Renormalization group analysis 

The Hamiltonian description in terms of the dilute Potts model allows a renormalization group 
analysis of percolative properties within the lattice gas. The case we are interested in is the 
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two-dimensional one, for which the magnetic and percolative transitions in the Ising model both 
take place at T c for H = 0. The renormalization group analysis for this case was performed in 
[U (see also \18\ 119]). We now recall the main results of this analysis casting them within the 
field theoretical language needed for our subsequent purposes. 

We look for fixed points of the Hamiltonian (J3|) in d = 2 for q — ► 1. Since the percolative 
properties do not affect the magnetic ones, we need to be at a magnetic fixed point to start 
with. The only non-trivial such fixed point is at T = T c , H = 0; there we look for fixed points 
of the residual coupling J. 

Before starting this search, let us recall that in two dimensions one can associate to a fixed 
point of the renormalization group a number c, called central charge, which grows with the 
number of dynamical degrees of freedom [201 EH I22j . A fixed point (i.e. a conformal field 
theory) with central charge c parameterized as 



together with infinitely many other less relevant fields (called descendants) whose dimensions 
exceed ()10p by integers. The field (pi^i with dimension is the identity /. For m integer and 
larger than 1 the theory admits a 'minimal' realization in which the operator product expansion 
closes on the primaries with r and s positive integers up to m — 1 and m, respectively, and 
their descendants [20j . Considering that within this set of primaries each dimension appears 
twice and excluding the identity, the number of independent and non-trivial primaries within a 
minimal model is m(m — l)/2 — 1. 

Back to our problem, since no dynamical degrees of freedom are associated to the percolative 
properties, for all the fixed points with T = T c , H = the central charge is that of the critical 
Ising model, i.e. c = 1/2 (m = 3). Clearly, the decoupling point J = yields a trivial fixed 
point at which percolation plays no role (ps = 0) and we deal with the pure magnetic fixed 
point of the Ising model. The latter is well known to be described by the minimal realization 
of the m = 3 conformal theory, which indeed contains two relevant scalar fields: the spin field 
a with dimension X a = X\^ = 1/8, and the thermal field e with dimension X £ = $ = 1. It 
follows that at this fixed point the coupling J is conjugated to an irrelevant field. 

In order to progress, we need to recall that the dilute Potts model (j4|) (with T which can 
be fixed to T c for the time being) admits two distinct lines of fixed points for q < 4 |23j. The 
first one is the critical line of the undilute (A = +oo) Potts model [24] . while the second is a 
tricritical line arising at some A c yielding a finite concentration of vacancies. Both lines are 
described by conformal theories with central charge ([9]) and the following relations between q 



C 



(9) 



c= 1 - 



m(m + 1) 



contains scalar fields (p r>s (called primaries) with scaling dimension 




(10) 
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and m [251 [26] 



y/q = 2 sin 



2(t + l) ' 



for the critical line . 



t + 1 for the tricritical line . 



(11) 



The two fixed lines meet at g = 4 (m = oo, c = 1) and the transition is first order above this 
value. The Potts spin field has scaling dimension X s (q) which coincides with -X~( m -i)/2,(m+i)/2 
along the critical line, and with ^ m /2, m /2 along the tricritical one. As for the fields invariant 
under the S q symmetry of color permutation, the leading scaling dimension Xt 1 (q) along the 
critical (tricritical) line is A^i (^1,2); since also the dilution is relevant at tricriticality, a second 
relevant, S^-invariant scaling field is present along the tricritical line and corresponds to ^1,3, 
with dimension X t2 {q) = -X"i,3- 

In the limit q — > 1 relevant for the percolative properties, the critical line gives m = 2, c = 0, 
in agreement with the fact that the undilute case corresponds to random percolation, which 
carries no dynamical degrees of freedom. Hence, the critical exponents for random percolation 
are determined by those of the q — ► 1 pure Potts model. In particular, the cluster size exponent 
coincides with the Potts susceptibility exponent 7. This and the correlation length exponent v 
are given b}{§ 

1/y = 2 - X u , 1 /v = 2- 2X S , (12) 

and at the random percolation point are determined by the values X s (l) = 5/48 and Xt 1 (l) = 
5/4 on the critical line. These values characterize the infinite temperature fixed point of Fig. 1. 
Notice that the minimal realization of the m = 2 fixed point contains only the identity: critical 
percolative properties are not described by minimal conformal field theories. 

The limit q — ► 1 along the tricritical line gives m = 3, c = 1/2. This is also expected because 
at q = 1 there are no Potts degrees of freedom and the critical degrees of freedom are those 
of the lattice gas, which is an Ising model. As a consequence, within the model (J2|) this fixed 
point corresponds to T = T c , H = and to some J = J*. Contrary to the case J = discussed 
above, it is a non-trivial fixed point for correlated percolation. The percolative exponents are 
determined here by the tricritical values X s (l) = 5/96, A" tl (l) = 1/8 and Xt 2 (l) = 1. Once 
again, the first of these dimensions does not belong to the minimal realization of the m = 3 
theory. 

Notice that the dimensions of the the two relevant S^-invariant fields on the tricritical line 
coincide at q = 1 with those of the Ising spin and thermal fields, which are conjugated to H and 
T, respectively. As a consequence, the remaining coupling J in is irrelevant at J*. Since 
it was irrelevant also at J = 0, consistence of the renormalization group flows requires a third, 
intermediate fixed point where J is relevant (Fig. 3). This was located in [4] at J = 2/T c as a 



2 We recall that the values of the critical exponents depend on the direction in coupling space along which 
the fixed point is approached. In (|12|l we take into account that on the tricritical line there are two 5 ? -invariant 
relevant fields. 
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Figure 3: Renormalization group flows in the coupling J for the Hamiltonian @ in two dimen- 
sions with q — > 1, T = T c , H = 0. While J = is just the Ising magnetic fixed point, 2/T c and 
J* are also percolation points for Ising droplets and Ising clusters, respectively. 

consequence of the identity 

-^lj= 2 /r = |E(^^- 1 ) + ( ln ^- 2 ^)E^o, ^ = 0,1,..., q. (13) 

This equation is easily checked using U = 8 ai ,+, 8 auCJ] = 8 ai ,+8 ajt+ + 8 ai ,-8 aj - and 8 ai -8 aj - + 
8si,sj8ai,+8aj,+ — &m,iJ,j-i where is a site variable which takes the value Si if <Tj = 1, and if 
CTj = —1. Substitution into the original form of the Hamiltonian gives f)13|) with instead of 
Ui and without the lng term. The latter arises because the q ways of coloring a site with spin 
down produce a factor q Nv when switching from the partition sum over {fii} to that over {vi). 
The Hamiltonian (| 13j) is that of a (g + l)-state Potts model and is critical for 2H = Inq and 
some g-dependent value of T. As q — > 1 we obtain an Ising model with critical point at H = 
and T = T c . 

It is clear from (|13p that J = 2/T is a special case of Q in which the lattice gas variable and 
the Potts spin are treated symmetrically, so that they have the same scaling dimension at the 
fixed point J = 2/T c : X s = X a = 1/8; similarly, = X £ = 1. This means that the KF clusters 
with pb = 1 — e~ 2 l T percolate at T c with exponents 7 and v which coincide with those of the 
magnetic susceptibility and correlation length in the Ising model. Since these are requirements 
of the droplet model [5] meant to describe magnetic transitions in a cluster language, the KF 
clusters with ps = 1 — e _2//T are called Coniglio-Klein Ising droplets [I]. 

Concerning Ising clusters, the renormalization group pattern of Fig. 3 finally leads to the 
conclusion that their critical exponents at T = T c , H = are those of the fixed point J*, onto 
which J = 00 (pb = 1) renormalizes at large distances [3]. The mean cluster size exponent 
in the limit T — > T~, H = was first evaluated by series expansions in [12], with the result 
7 ~ 1.91, quite close to the exact one 7 = 91/48 = 1.895.. coming from X s = 5/96 and v = 1. 
This exact value was identified in |19j . 

We finish this section recalling that at the percolation point clusters behave as fractals, i.e. 
their size grows with the linear extension L as L D , with a fractal dimension D smaller than the 
space dimensionality d [JJ. The fractal dimension is easily determined considering the number 
Soo(L) of sites belonging to the incipient infinite cluster inside a box of side L. This goes like 
PL d , where P is the density of the infinite cluster, namely the percolative order parameter, 
which goes like the Potts magnetization per site, i.e. like L~ Xs . Hence D = d — X s . In<i = 2 the 
values of X s given above determine the fractal dimensions 91/48 = 1.89.. for random percolation, 
187/96 = 1.94.. for Ising clusters, and 15/8 = 1.87.. for Ising droplets, which are less dense than 
Ising clusters, as expected. 
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Figure 4: Phase diagram of the field theory (|15p describing the scaling dilute Potts model. The 
ferromagnetic phase transition surface corresponds to A = 0. On this surface, the trajectories 
originating from the tricritical line are massive for g < (first order transition) ; they are massles 
and flow into the critical line for g > (second order transition). 



4 Field theory of the scaling limit 

The scaling region of lattice models around second order phase transition points can be described 
in a continuous, field theoretical framework. With the notation for the fields introduced in the 
previous section, the scaling limit of the Hamiltonian (JTJ) in d = 2 is described by the Ising field 
theory ([27] for a review) with action 



Aising = J^ c — t J d 2 xe(x) — h J d 2 xa(x) 



(14) 



where is the action of the conformal field theory with central charge c = 1/2, and 

r ~ T — T c , h ~ H as the critical point is approached. The action (I14p encodes all the 
universal features of the ferromagnetic phase transition. The universal percolative properties of 
the Ising model are instead contained in the scaling limits of the dilute Potts Hamiltonian 
with q — > 1. It follows from the discussion of the previous section that, as far as Ising clusters 
are concerned, the renormalization group trajectories which matter are those originating from 
the line of tricritical fixed points. These are described by the field theory 



a a trier 



g j d 2 xipi- i {x) - A / d 2 x(fi^{x) 



(15) 



where A^p^p is the action of the conformal theory with central charge ([9]), related to q by 



(|TT|) with t = m — 1, and, as we saw, ip± t 2 and ipip are the two relevant S^-invariant fields at 
tricriticality. 

Let us first discuss some general features of the field theory (I15p for 1 < q < 4. It follows 
from general results [28] that the theory is integrable when at least one of the two couplings g 
and A vanishes. Integrability allows to establish that the critical surface separating magnetically 
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Figure 5: Qualitative behavior (along the dotted path on the right) of the magnetic order 
parameter M q and of U q = (921,2) in the scaling dilute Potts model (|15|) . For q = 1 this model 
becomes the Ising field theory (|14p with g = r and A = h; M\ and XJ\ become the Ising 
percolative and magnetic order parameter, respectively. 



ordered and disordered regions of the scaling dilute Potts model corresponds to A = (Fig. 4): 
for A = 0, g < 0, (|15p describes the first order part of the transition surface on which the q 
ordered ground states are degenerate with the disordered one [29]; for A = 0, g > 0, instead, the 
theory is massless [2TJ [30] and describes the second order part of the transition surface, spanned 
by the trajectories flowing from the tricritical to the critical line of fixed points. 

For A < 0, the theory (I15h possesses a unique vacuum \Q) corresponding to the disordered 
ground state; for A > the degenerate vacua |fii), i = l,...,q, correspond to the q ordered 
ground states. The order parameter of the ferromagnetic transition is the Potts magnetization 
per site M q , i.e. the expectation value of the field S s (x),j ~ taken over for A < 0, and 
over |fL) for A > 0. At the first order transition the q + 1 vacua are degenerate, i.e. (971,3) is 
the same on all vacua. The expectation value U q = (^1,2), like M q , is discontinuous across the 
first order transition surface, its discontinuity corresponding to the latent heat [29]. Both M q 
and U q vanish on the second order transition surface (Fig. 5). 

We are now ready to discuss the limit we are actually interested in, q — ► 1 + , which we 
know a priori is peculiar. Indeed, it follows from the scaling dimensions given in theprevious 
section that, as expected, A\ = Aising, with g = r and A = h. Since Ai s i ng is massivqj for any 
value of the couplings away from r = h = 0, the massless trajectories of (|15p for A = 0, g > 
have to become massive as q — > 1. Remarkably, this phenomenon can be described analytically 
exploiting the integrability of (|15p with A = 0. 

Let us recall that use of (non-conformal) integrability in two-dimensional field theory requires 
switching to a particle language. Indeed, integrable field theories are solved by exact determina- 
tion of the ^-matrix of the associated scattering theory in (1 + l)-dimensional space-time [31] . 
In the ordered phase A > 0, the elementary excitations of (|15p on which the scattering theory is 

3 Field theoretically, due to the reflection positivity of (|14j> . this conclusion follows from Zamolodchikov's c- 
theorem |21) . 



11 



built are the q(q — l)/2 kinks which interpolate between the q degenerate vacua, together with 
their antikinks; due to S q symmetry all these kinks have the same mass |32j . In the disordered 
phase A < 0, instead, the elementary excitations are ordinary particles A^, k = 1, ... ,q — 1, 
forming a degenerate multiplet in which the antiparticle coincides with A q _k [33] . The mass 
of all these excitations goes to zero when the second order transition surface is approached 
(A — > 0, g > 0). Since the q vacua of the spontaneously broken phase coalesce in this limit, the 
elementary excitations on the transition surface are to be identified with the massless limit of 
the particles A^, which, in (1 + 1) dimensions, are right/left movers with energy p° = (/i/2)e =t6 ' 
and momentum p 1 = ±p°, fi being a mass scale and a rapidity parameter. 

While the number of these massless degrees of freedom clearly vanishes as q — > 1, the 
emergence of a massive particle in the same limit follows from the study of the right-left scattering 
amplitudes^. The latter are known exactly [30] within a particle basis which is not the Potts 
basis we are discussing. This is, however, immaterial as far as the analytic properties we are 
interested in are concerned (see e.g. [37]), and the only think we need to know is that the poles 
of the amplitudes are determined by the factor [30] 



1 

exp 



cosh p(iir — 6) 



"°° dx sinh § 9x 

' l I T ~~C~x — ± t-^B3H — 

/ x sinn cosh | tt 



(16) 



where p = l/(m — 1) and 9 determines the square of the center of mass energy as s = p 2 e 8 . 
As usual, the scattering amplitudes have a multi-sheet structure in the complex s-plane. The 
'physical sheet' on this plane corresponds to Im.9 G (0, 7r), and the 'second sheet' to Im# € 
(0, — tt). Writing the exponential part of (fT6j) as 5 , _ 1 / 2 (^)/5'i/2(^)) with 



QQ r(l + (2n + |- 7 ),-g)r(i + (2n + l- 7 )p + ^ 
iioT a + (2n + §- 7 ) P + ^ T (i + (2n+*- 7 )p-g 



one can see that for q > 1 (i.e. m > 3) (fT6j) has no poles on the physical sheet, and possesses 
a single pole, coming from the 1/ cosh prefactor, at #o = — ivr(m — 3)/2, which is on the second 
sheet for m G (3,5). This pole at sq = p 2 e e ° on the second sheet (Fig. 6) corresponds to a 
resonant particle A in the right-left scattering channel. When q — > 1 + , the pole approaches the 
physical sheet and gives a narrow resonance with square mass Re so — p 2 and inverse lifetime 
Imso//U (x (q — l)p, which coexists with the q — 1 massless particles. At q = 1, no massless 
particles are left and A provides the stable massive particle of the Ising field theory (I14p with 
h = 0, r > 0. 

Since the number of components of Sg-multiplets vanishes as q — > 1, the particle A we are 
talking about must be an ^-singlet, i.e. arises in the right-left scattering channels A^Ag-k- 
This prevents the correlators of S^-invariant fields from vanishing as q — > 1. Indeed, let us 
denote by 4> such a field (scalar, for simplicity) and consider the spectral decomposition of its 

4 See |34l 135] for the theory of integrable massless scattering. 

5 See [31] for a discussion on 5 q -invariant scattering theories and change of particle basis. 
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Figure 6: Location sq of the resonant pole on the second sheet of the s-plane. The pole describes 
an anticlockwise motion along the dashed path as m decreases from 5 to 3: sq is purely imaginary 
at m = 4 (q = 2) and real and positive at m = 3 (q = 1), where the associated particle becomes 
stable. 



two-point function G(x) = (4>(x) (f>(0)) over asymptotic particle states, in the limit q — ► 1. The 
contribution to G{x) with the lowest number of intermediate particles is 

lim(g-l) / de 1( %Kn|0(O)|A fc (ei)Vfc(^)>| 2 e- |!B|Ba ' o(9l ' fl3) , (18) 

where the prefactor q — 1 comes from the summation over k, and E2, 0(^1,^2) = /^(e 6 * 1 + e~ e ' 2 )/2 
is the energy of a right-left pair. Right-right and left-left pairs do not contribute to the limit 
because only the right-left matrix element inherits from the right-left scattering amplitude the 
resonance pole at 0\ — 62 = 80, with #0 oc —i(q — 1) for q — > 1. Hence, calling the residue of 
the matrix element on this pole, (|18[) can be written aqj 

lim (q-1) [ d(3d0 - j/^' 2 - e -Wifc.o((/J+9)/2,(/?-*)/2) K f df3 e -\x\E h »Q3) ; (19) 
g->i 7 (0 — 9q)(8 + 0o) J 

where E\^{(3) = /zcosh(/3/2) is the energy of a particle with mass fj,. We see in this way 
how correlations mediated by massless particles on the second order surface are replaced by 
correlations mediated by the particle A along the massive trajectory at q = 1 (Fig. 7). Fields 
which are not S^-invariant do not couple to A and have zero correlations (i.e. are absent) at 
q = 1 (as required by S\ = I). Recalling that the ^-invariant fields along the tricritical Potts 
line are /, 921,2, ¥>i,3 an d their descendants, which become /, a, e and their descendants at q = 1, 
the mechanism we illustrated explains in particular how the theory becomes minimal for this 
value of q. 

Summarizing, the percolation line Hq{T) of Fig. 1 is mapped in the scaling limit onto the 
case h = of the field theory (fl4"]l. More precisely, as expected from ([7J), this field theory 
describes only the magnetic properties of the Ising model; the universal percolative properties 
of Ising clusters are described by the 'embedding' of (|14p into the dilute Potts field theory (|15p 
with q — > 1. In this limit, the expectation values M q and U q shown in Fig. 5 become the Ising 



By relativistic invariance, the two-particle matrix element of a scalar field depends on the rapidity difference, 
and the one-particle matrix element is rapidity-independent. 
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Figure 7: Correlations mediated by two massless particles at q > 1 are mediated by one massive 
particle at q = 1. Fields (j> which are not S^-invariant do not couple to the massive particle and 
are absent at q = 1. 

percolative and magnetic order parameter, respectively. While for r < both transition are 
first order, for r > there is only a continuous percolative transition, a circumstance which is 
explained analytically by the evolution of the resonance pole discussed above. 

Notice that for the transition at r < the situation is simpler. On the first order part of 
the transition surface of (|15p the q ordered vacua |f2j) are degenerate with the disordered one 
| Q). The elementary excitations are q kinks interpolating from \Q) to together with their 
antikinks [29] . When q — > 1 one recovers straightforwardly the single kink of the Ising model. 



5 Universal scaling limit of the Kertesz line 

We saw in section 3 that for H = the Coniglio-Klein droplets, i.e. the KF clusters with 
Pb = 1 — e~ 2 / T , percolate at T c (also in d = 3) with critical exponents coinciding with the 
magnetic exponents. Also, contrary to the case of Ising clusters, an infinite droplet above T c in 
zero field is excluded by the vanishing of pb at infinite termperature. It is easy to see, however, 
that when H ^ the critical properties of the magnetic and percolative degrees of freedom 
no longer coincide. Indeed, for H = +oo all the spins are up and we are left with a random 
bond percolation problem with occupation probability ps, which is critical for some value p B . 
This means that there is a line Tk(H), called the Kertesz line [6], going from Tr-(O) = T c to 
Tfc-(+oo) = — 2/ln(l —p B ), which is a percolation line for the Coniglio-Klein droplets (Fig. 1). 

The universal features of these droplets are described by the scaling limit of the Hamiltonian 
(|13p . namely, in d = 2, by the field theory 

■Adroplets = -^CFT ~ T Q J ^ X ^2,1 0*0 + 2/lg / d 2 X 5 u{x)fi , v(x) = 0, 1, . . . , q , (20) 

where Agp^ accounts for the critical line of a (q+l)-state Potts model, r q measures the deviation 
from the critical Potts temperature, and h q is a magnetic field pointing in the u = direction. 
Since in two dimensions the Potts transition is second order as long as the number of states does 
not exceed 4, the above action is intended for q < 3. For q = 1, (|20|) with t\ = r and h± = h 
gives back the Ising field theory (|14p . 

The renormalization group trajectories flowing out of the fixed point at r q = h q = for 
h q > are labelled by the dimensionless parameter 

/, (2— Jf2,i)/(2— X a ) / n1 \ 
T]q = Tqjh\ (21) 

where the thermal and magnetic scaling dimensions X2 t i and X s are those given in section 3 
for the Potts critical line, up to the substitution q — > q + 1; rji = rj = r/h 8 / 15 labels the Ising 
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Figure 8: Scaling limit of Fig. 1 within the parameter space of the Ising field theory (j!4|) . The 
percolative order parameter is non-zero in the whole upper half-plane for Ising clusters, and to 
the left of the Kertesz trajectory tjk for Coniglio-Klein Ising droplets. 



trajectories. For h q = +oo the state v = is forbidden and one is left with a q-state Potts model. 
This means that for 1 < q < 3 there is a critical trajectory rf q which flows from the (q + 1)- 
to the q-state Potts fixed point. The S^-invariance of (f20j) is spontaneously broken for tl, < rjg. 
Since the Potts Curie temperature increases as q decreases, r(i is expected to be positivqll 

The massless trajectories r\ c q span, as a function of q, a critical surface which plays for the 
droplets exactly the same role the second order surface considered in the previous section played 
for the Ising clusters. In the present case, however, while it can be seen that the number of 
massless particles indeed vanishes as q — > 1 [33], lack of integrability does not allow to follow 
analytically the origin of the mass gap at q = 1. The most likely mechanism remains the one seen 
in the previous section, namely at least one neutral massive resonance which becomes stable at 
(7 = 1. The possibility that this neutral particle is stable already for q > 1 cannot be excluded, 
but requires a mechanism which prevents the decay into the massless excitations. 

Either way, it follows from the same line of arguments developed for the clusters that the 
value T}\ = t)k of the Ising parameter r/ determines the universal scaling limit of the Kertesz 
line (Fig. 8). An infinite Coniglio-Klein droplet exists in the sector of the r-h plane with h > 0, 
r] < t]k- Percolative critical exponents measured along the trajectory t\k (as along the trajectory 
r\ = +oo for clusters) are those of the percolative infrared fixed point, i.e. the random percolation 
fixed point with central charge c = 0. 

Numerical data for the Kertesz line of the square lattice Ising model in the vicinity of the 
magnetic critical point are given in [38]. Using them, together with the known relations between 
lattice and continuum parameters (see e.g. [39]), we find_| rjx — 0.12. We recall that, as a result 
of a series of theoretical and numerical studies (see [27] for references), very much is known 

7 The field theory (|20[) has been considered in [33] with the purpose of describing the qualitative evolution of 
the particle spectrum in parameter space. The actual value of rf q was not essential there and was naively identified 
with zero. 

8 This and the subsequent values of 77 refer to the normalization 

lim |*| 1/4 <a(*)a(0)) = lim \x\ 2 (e(x)s(O)) = 1 (22) 

of the fields in CH]). 
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about the 77-dependence of the particle spectrum of the field theory (TH| . In particular, the 
theory possesses a single stable particle for 77 > rjm , and two stable particles for 77/2) > f] > T)(3) ; 
with r?(2) — 0.33 and 77(3) ~ 0.022 [ID] . It follows that r]x falls inside the second region, so that, 
if the resonance scenario for the production of the mass gap as q — > 1 along the critical surface 
applies also to this case, two particles need to become simultaneously stable at q = 1. 

6 Conclusion 

In this paper we showed explicitly the non-perturbative analytic mechanism through which 
clusters of up spins in the two-dimensional Ising model undergo a second order percolation 
transition at values of temperature T and magnetic field H for which there are no thermodynamic 
singularities. The Kasteleyn-Fortuin formulation of percolation leads to consider a dilute (/-state 
Potts model in which the Ising degrees of freedom are associated to the dilution and the Potts 
spins are auxiliary site variables. In the three-dimensional parameter space of T, H and q there 
is for q > 1 a surface of second order phase transition which, in the scaling limit, corresponds to 
an integrable field theory. This allows to show that the q — 1 massless particles on the surface 
produce in their scattering a massive resonance whose lifetime is proportional to l/{q— 1). When 
the limit q — > 1, needed to get rid of the auxiliary Potts variables, is taken, the percolative 
correlations, which are determined at first order in q— 1, are mediated by the massless particles, 
giving rise to the second order transition line. The thermodynamic observables, instead, are 
determined at q strictly equal to 1, where no massless particle is left and the correlations are 
mediated by the massive particle which becomes stable. The fact that the resonant particle 
is a singlet under the Potts permutational symmetry S q leads to a smooth transition from the 
massless to the massive regime for the thermodynamic observables as q — ► 1. 

We have analyzed in the same framework also the Coniglio-Klein droplets, which are obtained 
from the Ising clusters by a partial, temperature-dependent depletion and produce at the Curie 
point percolative exponents which coincide with the magnetic ones. The corresponding field 
theory in the Kasteleyn-Fortuin representation describes the scaling limit of a (q + l)-state 
Potts model with a magnetic field which reduces the symmetry to S q . Again, for q > 1 there is 
a massless surface bounded at q = 1 by a massive trajectory which corresponds to the universal 
scaling limit of the second order percolation line for the droplets (the Kertesz line). In this 
case the massless surface is not integrable, but a mechanism analogous to that observed for the 
clusters is likely to account for mass generation at q = 1. 

On the lattice, the second order percolative line for the clusters (droplets) goes from the Curie 
point to a finite value of the magnetic field (temperature) at infinite temperature (magnetic 
field). The asymptotic value is determined by the critical probability of random site (bond) 
percolation on the given lattice, and is non-universal. Our field theoretical description of the 
scaling limit retains only the universal features. The second order critical lines correspond 
to renormalization group trajectories flowing from the correlated percolation fixed point, with 
central charge c = 1/2 and different dimensions of the percolative order parameter for clusters 
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and droplets, to the random percolation fixed point with c = 0. These universal trajectories are 
identified by two values of the Ising field theory parameter r/ = r//i 8 / 15 : rj = +00 for the clusters 
and rj ~ 0.12 for the droplets. The latter value is obtained from existing numerical data and 
implies that the Kertesz trajectory falls in the sector of the Ising field theory with two stable 
massive particles. 
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